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Abstract 

In the continuum, the single flavor massless Schwinger model has an exact global axial 
U{1) symmetry in the sector of perturbative gauge fields. This symmetry is explicitly 
broken by gauge fields with nonzero topological charge inducing a nonzero expectation 
value for the bilinear i/jiIj. We show that a lattice formulation of this model, using the 
overlap formalism to treat the massless fermions, explicitly exhibits this phenomenon. A 
Monte Carlo simulation of the complete system yields the correct value of the fermion 
condensate and shows unambiguously that it originates from the sector of topological 
charge equal to unity. 
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Introduction 



In chiral gauge theories fermions appear bilinearly in the action and interact with each 
other only by the exchange of gauge bosons. The bihnearity in fermion fields is preserved 
in most methods of rcgularization. Then the regularization becomes defined for any single 
(complex) irreducible representation carried by the chiral fermions, be it anomalous or 
not. Using this definition one can construct, in particular, vector gauge theories where the 
multiplets are assembled in left-handed/right-handed pairs per representation. Massless 
QCD is an example of such a vector gauge theory and it has an anomalous global U{1) 
symmetry that is violated by topology induced 't Hooft vertices [1]. Therefore, any reg- 
ularization of chiral fermions that preserves bihnearity should break the anomalous U{1), 
but conserve the nonanomalous global symmetries. Some of the existing proposals may 
have problems in this respect; we know of no examples where this point has been tested 
with quantum fluctuations of gauge bosons taken into account. 

Here we subject the overlap formalism to this test. For smooth, classical, gauge 
fields we already know that the overlap formalism works correctly [2,3]. The purpose of 
this paper is to report a Monte Carlo simulation of the single flavor, strictly massless 
Schwinger model [4] formulated on the lattice. Each of the two chiral components of the 
Dirac fermion are realized using the overlap formula. Of course, this model can also be 
formulated on the lattice using conventional techniques [5] , but then a massless limit has 
to be taken, there is difficulty in realizing a single flavor, there is no factorization of the 
theory into left and right sectors for flxed gauge flelds, and the role of topology is obscured. 
Anyhow, our main point is not to present yet another solution of the Schwinger model, 
but, to establish a useful test for proposals to regularize chiral gauge theories and show 
that the overlap passes this test successfully. 

In order to make it absolutely clear that the overlap succeeds in computing something 
that no previous approach to regularizing chiral fermions could, we picked as our observable 
just the 't Hooft vertex which together with its hermitian conjugate makes up the 

mass bilinear ■^'^ more familiar in this context [6]. Our observable receives contributions 
only from gauge flelds of topological charge one. The Schwinger model is a good model to 
pick because it has the necessary physics ingredients, its study is not particular demanding 
computationally and, above all, it is exactly soluble in the continuum, even for flnite size 
systems. Various aspects of a two flavor version of the model have been investigated in 
detail in recent years both in the continuum and on the lattice by Azakov, Joos and Dilger 
[7] . Here we rely mostly on the continuum analysis of Sachs and Wipf [6] . 
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The model 



The continuum action in Euchdean space is: 

S'[V5,V',-4^] = {(^tiud^Auf - j i>Lcr^^{d^+iA^)ij;L- J i'R(J*^{df^ + iA^)'ij;R 

where /U, z/ = 1, 2, ai = 1 and o"2 = i- iI'r,Lti1^r,l are Grassmann variables and is a 
U{1) gauge field. The model is defined on a torus of fixed physical size / x / with periodic 
boundary conditions for the electric field, E = e^^d^A^, and the fermions. The topological 
charge q is an integer given hj q — J E. The expectation value of the fermion bilinear 
is 



The lattice formulation begins by replacing the continuum torus with a toroidal square 
lattice consisting of L sites in both directions. The lattice spacing a is a = l/L. Gauge 
fields U^{x) = e'°"^'^(^) are associated with the links of the lattice. The overlap formalism 
[2] is used to deal with the chiral fermions on the lattice. This results in the following 
expressions for N and D appearing in the expectation value of the fermion bilinear: 



N 



= ^2 I] / n dUf,{x)e-^^^^^ [\u<L- \a^a\L+ >u\^ + \u<L- \aljL+ >u \^] 

xa 

D = jl[dU^{x)e-^^^'''^\u<L-\L+ >u\'' 

The a-term in N comes from gauge fields carrying topological charge +1 and is the overlap 
expression for < iPri/jl > while the a^-term comes from charge —1 and represents < 
'4^L'<pR >• In our simulations we could have looked at only one of these terms; we took 
both in order to increase statistics. Analytically, the two overlap matrix elements have 
exactly equal magnitudes. The CP violating angle 6 is set to zero and therefore only the 
square of the absolute value appears in the above formula. The states |L± >[/ are the 
ground states of many body Hamiltonians 

xa,y(3 

a, P = 1,2 and x = {xi,X2) with x^ = 0,1,2, ,L — 1. The single particle hermitian 
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hamiltonians H"^ are given by: 




The parameter m can be chosen in the range < m < 1 and is set to m = 0.9 here. Note 
that topological charge ±1 means that the number of single particle states filled in \L+ >u 
and \L— >u differ by ±1. 

We need a gauge action and an associated Monte Carlo method that will give reason- 
able estimates for the topological charge distribution. For this reason, we would like to 
avoid correlations between configurations as far as possible. With this in mind, we pick a 
single plaquette action. 



for the gauge action appearing in the formulae for and D. Up is the oriented product 
of the four link elements making up the plaquette p. For such an action the plaquettes 
are almost independent variables in the pure gauge theory. To minimize the finite lattice 
spacing effects arising from the gauge action, we chose the "heat- kernel" action [8]: 



This choice assures exact scaling of the string tension with the coupling g in the pure gauge 
theory. 



To compute the bilinear, we use a Monte Carlo method where configurations are drawn 
independently of each other. This avoids correlations between sequential samplings. To 



p 




The Monte Carlo method 



this end we replace N by N' and D by D' with N' 



= f and = § 



where 
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Gauge configurations will be generated stochastically using the distribution e ^3^^^ and, 
denoting averages w.r.t. e~^9^^^ by << .. >>, we evaluate probabilistically 

N' = J2J2« e-^«(^)+^^(^) [\u< L - |a^,a|L+ >u |' + \u< L - \aljL+ >u » 

xa 

D' =« e-^«(^)+^^(^)|t;< L- |L+ >u\^ » . 

Sg{U) is picked to be 

S'giU) = S,iU) - s,iU,*) 

where p* is some arbitrarily chosen plaquette. The probability distribution e'^f*^^-' factorizes 
into one factor of e"''''*-^*'-' for each p ^ p* . To generate a configuration, one draws a set of 
(L^ — 1) Up variables, each according to the plaquette factor e~^3^^p\ This fixes (L^ — 1) 
out of the (L^ + 1) link variables remaining after choosing a gauge on the L x L lattice. 
Two more link variables determining Polyakov loops winding around the torus are drawn 
randomly since the gauge action is unaffected by these loops. Since we are computing 
gauge invariant quantities the exact nature of gauge fixing is irrelevant and there is no 
need to average over gauge transformations. 



The parameters 



Having described the Monte Carlo method, we now proceed to determine the param- 
eters of the model. Since we are treating the overlap factors as observables in a pure gauge 
theory, this method will work only if the fluctuations in these factors are not too rapid. 
Technically, they are of the order of e*^^^ and since we are simulating a system of finite 
physical size I, c is also dependent on L. Based on the exact solution in the continuum we 
choose I so that the fluctuations in the fermionic determinant are small. We feel confldent 
in using the continuum results for our estimates since we have already shown that the 
overlap agrees with the continuum very well for sufliciently smooth flelds [2,3]. The zero 
momentum mode of the electrical fleld is quantized by topology, so the dependence on it 
is special. The fermionic determinant also depends on the zero momentum component of 
the gauge flelds, but these are only two variables and the dependence though sizable is 
not too violent. Previous work assures us that the overlap will follow continuum behavior 
faithfully also for these variables [9]. There remain now the nonzero momentum modes of 
the electrical fleld E{k) = (j){k)k'^. The bosonic action per mode is e i^''^^'^^' ^ yielding 
The fermion determinant factorizes among the modes [2] and, per mode, is 



5 



2 2 

given by e~2^Fl<^('^)l ^ ki e~k^ where /x = The smallest nonzero momentum is = ^ 

so we get a typical contribution of the order e . For /i/ ?s 3, the typical number in 
the exponent is expected to be quite small. If we sum over the higher modes the result 
changes a little but we do not need to complete this exercise for our simple estimates. 

Apart from the overlap factor, we also have the factor e"'^''*-'^-*''"'^^*-'^-* = e"'^^*^'^^*^ 
with t/*, = Hp/p* ^p- should note that our method does not implement importance 
sampling for this factor. We were willing to accept this drawback in exchange for the 
unhindered generation of the topological charge distribution. Since our physical volume is 
fixed, the coupling becomes weaker as the lattice gets finer and e^^iUp) gg^g niore peaked. 
Its sampling however has a roughly constant and quite broad width, so, unlike in the 
usual case where statistical errors per lattice decrease with size due to self averaging, here 
the statistical error will stay more or less constant as a function of size and therefore the 
computing cost per given accuracy will increase dramatically with size. Although this can 
be remedied to some extent, it is not necessary for the computation presented here. 

Gauge configurations contributing to A^' have unit lattice topological charge and those 
that contribute to D' have zero lattice topological charge. In all other cases the respec- 
tive overlap vanishes identically. There is a probability p{0) = p to generate gauge field 
configurations of zero lattice topological charge. In a sample of n lattices there will be 
j ~ pn lattices contributing to D' . ^ is Poisson distributed with standard deviation 



cr(0) = y ■ In the continuum one has p{k) — Me . For /i/ ~ 3 for pure gauge 

(ignoring the single plaquette p* in the continuum) the decay at higher /c's is suflficiently 
rapid to ignore all but k = 0, ±1 as a first approximation for the expected errors. Hence 
there is a factor of -r^ in the ratio ^ • This factor will be distributed with a relative error 

l—p u 

5 = , To minimize it we would like to work with p ^ h. This is indeed the case 

for fil fa 3. 

In view of the above considerations, we picked fil = 3. The continuum limit is ap- 
proached by scaling the lattice coupling g with the lattice spacing via g = Simulations 
were performed on lattices ranging from L = 3 to L = 10. Since 5 !v ^ is a lower bound 
on the relative error in we see that in order to reach an accuracy of a few percent, we 
need a few thousand lattices at least. This is quite acceptable if the lattice sizes L can be 
kept small. 

The continuum result for the dimensionless quantity = f(fil) is known exactly 

from the work of Sachs and Wipf [6]: /(3) ^ 0.236 and /(oo) ?s .284. While we have no 
reason to try to extrapolate numerically to infinite physical size, we like the fact the we 
are working on a system that is reasonably large in physical terms. The finite size effect 
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is seen from the exact solution to be of the order of 20% and it is nice that we don't have 
to ehminate it numerically. 



Algorithm 

We used very direct and unsophisticated algorithms. Double precision was used 
throughout the algorithm. The diagonalization of H='= was done using Jacobi's method 
and the needed determinants evaluating the overlaps were computed using LU decom- 
position. Since half of the configurations require the evaluation of determinants the 
algorithm goes, asymptotically, as L^. In practice, the larger prefactor in front of the 
corresponding to Jacobi diagonalization versus the one corresponding to determinant 
computation seems to be important for our sizes and the algorithm scaled slightly worse 
than L^. The worst case, L = 10, took an average of 5 minutes CPU time per lattice on 
an 100 MHz a-workstation. No attempts were made to speed up the code. The results of 
this paper were obtained using the equivalent of a few hundred hours on such a machine. 



Results 

Our results are summarized in the table below where the numerical results are listed 
as a function of L. Errors were estimated in several ways and all gave consistent answers. 
We see that within statistical errors the correct answer is obtained already at L = 5, so 
most of the computer time we used was not really needed. We felt, however, that it was 
necessary to ascertain that nothing bizarre was happening when one went closer to the 
continuum limit. 



L 


/l(3) 


# of confs. 


3 


0.0812(32) 


5000 


4 


0.1638(66) 


5000 


5 


0.224(11) 


5000 


6 


0.223(13) 


5000 


7 


0.235(14) 


6000 


8 


0.238(12) 


8173 


10 


0.224(18) 


4919 


oo 


0.236 


Continuum 



Table: /l(a*0 = {^t)D various lattice spacings aX. jil = 3. 
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Conclusions and Outlook 



In this paper we liave sliown tliat tlie overlap formalism for chiral fermions passed 
a certain "vector test". This test can be used as a check for any other proposal that 
regularizes a chiral gauge theory. It is simple, not very costly to apply and, at present, we 
believe it would be a good selection method. It is hard for us to see how a bilinear method 
for the regularization of chiral gauge theories could succeed without passing this vector 
test. To be sure, let us stress that the numerical result we have obtained here should be 
obtainable by traditional methods [5] . In these methods one does not have exact left right 
factorization. When one starts from a chiral method this factorization is built in. 

The next step would be to study a chiral Schwinger model in the overlap formalism. 
Here we do not have exact gauge invariance on the lattice in the anomaly free case although 
this is achieved in the continuum for smooth gauge fields. We have to check whether the 
Foerster-Nielsen-Ninomya [10] mechanism of gauge invariance restoration in the contin- 
uum limit indeed works on the lattice as we conjectured in [2]. Preliminary studies we 
have done are encouraging. Then the road to exhibiting fermion number violation on the 
lattice is open. 
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